THE ORTHOGONAL WEINGARTEN FORMULA IN COMPACT 

FORM 

TEODOR BANICA 

Abstract. We present a compact formulation of the orthogonal Weingarten formula, 
with the traditional quantity . . . ,i2k ■ ji, ■ ■ ■ ,j2k) — Jq^ hk replaced 

by the more advanced quantity I{a) = nit°-^ du, depending on a matrix of exponents 
a G M„(N). Among consequences, we establish a number of basic facts regarding the 
integrals /(a): vanishing conditions, possible poles, asymptotic behavior. 



Introduction 

The computation of polynomial integrals over the orthogonal group 0„ is a key problem 
in mathematical physics. These integrals are indeed known to appear in a wealth of 
concrete situations, coming from random matrices, lattice models, combinatorics. 

The standard approach to the computation of such integrals is via the Weingarten 
formula. This formula, originating from Weingarten's paper [12], and worked out by 
Collins in [B], then by Collins and Sniady in fS], is an identity of the following type: 

/ Wiiii • • • Ui^khk du = ^ 5^{i)5a{3)Wkn{T^, a) 

Here the sum is over all pairings of {1, ... ,2k}, also called Brauer diagrams [5], and 
the delta symbols, describing the coupling between indices and diagrams, are or 1. As 
for Wkn, this is the key combinatorial ingredient: the Weingarten function. 

The exact or asymptotic computation of Wkn is a quite subtle problem, and several 
results have been recently obtained on this subject. Let us mention here the work of 
Collins and Matsumoto [7], Matsumoto and Novak [TT] and Zinn- Justin [13], providing 
a deep insight into the combinatorics of Wkn- Also, the foundational part of theory has 
benefited from several abstract versions and generalizations, developed in [T], [H], jlUj . 

The starting point for the considerations in the present paper is the following elementary 
reformulation of the Weingarten formula: 



On i-l j^l 
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Here a G M„(N) is the matrix of exponents appearing in the original Weingarten 
formula, and the delta symbols are once again or 1. This formula is of course fully 
equivalent to the original Weingarten one, but the slight difference in the formulation 
leads to some potentially interesting consequences, that we will explore in this paper. 

The idea is that our matrix formulation, with a G M„(N) replacing the multi- index 
(ii, . . . ,i2k '■ ji, ■ ■ ■ ,j2k) € {1, • • • , n}^'', allows us to escape a bit from the heavy combina- 
torial context of [12], [6], |8], [7], [11], [13]. As an example, in the case where a has only 
a nonzero row, with even entries, the exact formula for our integral is as follows: 



This kind of result, while very basic, is barely visible with the original Weingarten 
formula, and provides an excellent motivation for our proposed a G M„(N) reformulation. 
Some other exact formulae of this type are actually available, see the comments below. 

So, let us go back now to the modified Weingarten formula, with the matrix a G M„(N) 
replacing the multi- indices This provides a combinatorial expression for the following 
integral, to be regarded from now on as being the "basic quantity" to be computed: 



In this paper we will browse through the various combinatorial formulae in [12], [6], 
[8], [Z], [H], [9], by converting them into concrete results about /(a). We will obtain 
in particular exact results regarding the vanishing, sign and poles of I (a), viewed as a 
rational function of n G N, and also about its asymptotic behavior, with n oo. 

We believe that all these results can be of help in connection with a key problem, 
namely the exact computation of I (a). This question, recently raised in [2] in connection 
with the Hadamard conjecture, obviously stands above the Weingarten technology, and 
of the state-of-art of group integral computation in general. However, we do believe that 
the problem is actually quite tractable, and a "magic formula" for /(a), in the case where 
a is zero outside its upper left 2x2 corner, is in preparation [3J. There is of course an 
enormous work waiting to be done, for passing from the 2x2 case to the n x n case, and 
it is our hope that the present results about I{a) will be of help in this process. 

The paper is organized as follows: 1-2 are preliminary sections, in 3-5 we state and 
prove the various Weingarten function results about I (a), and 6-7 contain some results 
on the possible poles of /(a), along with some concluding remarks. 
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Segala, J. Novak, J.-M. Schlenker and P. Zinn- Justin for various discussions on the subject. 
The author was supported by the ANR grants "Galoisint" and "Granma" . 





n 



n 



i=l j=l 



the orthogonal weingarten formula in compact form 3 

1. General considerations 

We are interested in the computation of arbitary polynomial integrals over the orthog- 
onal group On- These are best introduced in a "rectangular way" , as follows. 

Definition 1.1. Associated to any matrix a G Mpxg(N) is the integral 

with respect to the Haar measure of On, where n > p,q. 

We can of course complete our matrix with values, as to always deal with square 
matrices, a G Mri(N). However, the parameters p, q are very useful, because they measure 
the "complexity" of the problem, as shown for instance by the result below. 

Let x!! = {x — l){x — 3)(a; — 5) . . ., with the product ending at 1 or 2. 

Theorem 1.2. At p = 1 we have the formula 

. . (n-l)!!ai!!...a,!! 

'("^ ••■ (n + Ea.-l)!! 

where e = 1 if all are even, and e = if not. 

Proof. This follows from the fact that the first slice of 0„ is isomorphic to the real sphere 
^"-^ Indeed, this gives the following formula: 



/ (ai ... Og) = . . . Xg" dx 



This latter integral can be computed by using polar coordinates, and we obtain the 
formula in the statement. See e.g. [2j. □ 

Another instructive computation, as well of trigonometric nature, is the one at n = 2. 
We have here the following result, which completely solves the problem in this case. 

Theorem 1.3. At n = 2 we have the formula 

'a b\ {a + d)\\{b + c)\\ 



I 



c dj (^a + b + c + d+l)\\ 

where e = 1 if a,b,c,d are even, e = —1 is a,b,c,d are odd, and e = if not. 

Proof. When computing the integral over O2, we can restrict the integration to SO2 = S^, 
then further restrict the integration to the first quadrant. We get: 



c d 



TC 



[ {cos ty+'^ {sin tf+^dt 
Jo 



This gives the formula in the statement. □ 
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The above computations might tend to suggest that I (a) always decomposes as a prod- 
uct of factorials. This is far from being true, but in the 2x2 case it will be shown in ^ 
that /(a) decomposes as a (quite reasonable) sum of products of factorials. 



In this reminder of this paper we discuss the representation theory approach to the 
computation of I (a). The technology available here is quite advanced, and will lead to 
a wealth of concrete results about I (a), regarding its vanishing, theoretical poles, and 
asymptotic behaviour, and notably with a concrete formula for its first order term. 

Our main tool is a combinatorial formula for the polynomial integrals over 0„, whose 
origins go back to Weingarten's paper [12]. The treatment given here follows the papers 
|Bj, [8j, [7j, where this formula was fully formalized, and its true power revealed. 

The first remark is that the integrals over On of arbitrary polynomial quantities of type 
Mj^j^ . . . Ui^j^ vanish, unless s is even. In what follows, we focus on the s = 2k case. 

Let us recall that the pseudo-inverse of a real symmetric matrix G is the unique matrix 
W satisfying WGW = W and GWG = G. If G is invertible, then W = G-\ 

Theorem 2.1. We have the Weingarten formula 



where the various objects on the right are as follows: 

(1) The sum is over all pairings vr, a of the set {!,..., 2k}. 

(2) The delta symbols are or 1, depending on whether the indices i = {ii, . . . ,i2k) 
and j = (ji, . . . ,j2k) fit or not inside the corresponding pairings. 

(3) Wkn is the pseudo-inverse of the matrix Gkn{'^,o') = n'-"°P^^'^''^\ where we denote 
by loops{7r, a) the number of loops obtained when superposing vr, cr. 

Proof. Consider the matrix P G M„2fc(R) formed by all the integrals on the left, with 
k fixed and i,j varying. It follows from the general theory that P is the orthogonal 
projection onto the space Fix{u®'^^) of fixed vectors of the 2/c-th tensor power of u. 
By a well-known result of Brauer [5] the space Fix{u®'^^) is spanned by the vectors 
= Sj5^(i)eij (g) . . . (g) Cj^j,, with 71 ranging over all the pairings of {1, ... , 2k}. Now since 
the Gram matrix of these vectors is < ^n,C,a >= Gkn{T!',o'), when computing P we have 
to invert this Gram matrix, and we obtain the formula in the statement. □ 

We refer to [6], [8], [7] for full details regarding the above proof. Let us also mention 
that the Gram matrix Gkn is actually invertible for n > k. See [8]. 

Our first task is to convert the above formula, by using our compact notation I (a) for 
the polynomial integrals over On- We restrict attention to the case where Eajj is even. 



2. The Weingarten formula 
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Theorem 2.2. We have the Weingarten formula 

Ko) = 5^{ai)5a{ar)WkniT^, o") 

7r,(j 

where k = Sajj/2, and where the multi-indices ai/a^ are defined as follows: 

(1) Start with a G Mpxg(N), and replace each ij-entry by aij copies of i/j. 

(2) Read this matrix in the usual way, as to get the multi-indices ai/ar. 

Proof. This is simply a reformulation of Theorem 2.1. Indeed, according to our definitions, 
the integral in the statement is given by: 

/(a) = J Uu . . ^.Uu Ui2 . . jUi2^ ... Upq. . .Upq du 

an 112 apq 

Thus what we have here is an integral as in Theorem 2.1, the multi-indices being: 

ai = (Ij^^^^ v_;.,>ii ... p ... p) 




ar = ... g ... g) 

an ai2 

The result follows now from the Weingarten formula. □ 

We are now in position of deriving a first general corollary from our study. This extends 
the vanishing results appearing in Theorem 1.2 and Theorem 1.3. 

Corollary 2.3. We have I{a) = 0, unless the matrix a is "admissible" , in the sense that 
all p + q sums on its rows and columns are even numbers. 

Proof. Observe first that the left multi-index associated to a consists of fci = Soij copies 
of 1, A;2 = Sa2j copies of 2, and so on, up to kp = T,apj copies of p. In the case where one 
of these numbers is odd we have STr{a) = for any vr, and this gives /(a) = 0. 

A similar argument with the right multi-index associated to a shows that the sums on 
the columns of a must be even as well, and we are done. □ 



3. Basic asymptotic study 

A natural question now is whether the converse of Corollary 2.3 holds, and if so, the 
question of computing the sign of I{a) appears as well. These are both quite subtle 
questions, and we begin our investigations with a n oo study. 

The basic result here, known since [S], [H], states that the Weingarten matrix is asymp- 
totically diagonal, in the sense that we have Wkn{T^, cr) = n~^{5T,a + 0(n~-^/^)). 

We present below a complete proof for this fact, by following a slighly improved method, 
which gives a better estimate, along with some concrete bounds. 
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Theorem 3.1. The Weingarten matrix is asymptotically diagonal, in the sense that: 

WkniTT, a) = n-\6^^ + 0{n-^)) 
Moreover, the 0{n~^) reminder is asymptotically smaller that {2k/ eYn~^ . 

Proof. It is convenient, for the purposes of this proof, to drop the indices k,n. We know 
that the Gram matrix is given by G(7r, a) = n'-°"^^'^'^''^\ so we have: 



G(7r,a) 



for t: — a 
n, V?, . . . , n^~^ for TT 7^ (T 



Thus the Gram matrix is of the form G — n^{I -\- H), with ||i/||oo < ^- Now recall 
that for any K x K matrix X, we have the following lineup of standard inequalities: 

Halloo <||^||< 11^112 <i^||^|U 

In the case of our matrix H, the size is K = (2A;)!!, so we have \\H\\ < Kn~^. Now 
from {I + H)-^^I-H + H'^-H^ + . ..wegei 1 17 - (7 + | < ||/7||/(1 - ||/7||), so: 

Wl-nJ'WW^ = ||/-(l + i7)-i|U 

< ||/- + 

< 11^11/(1-11^11) 

< Kn-^/{l- Kn-^) 

= K/{n-K) 

Together with the standard estimate K k, {2k/ e)^, this gives the result. □ 
Theorem 3.2. We have the estimate 

I{a)^n-' (flfla.j\\ + 0{n-'] 

\i=l 0=1 

when all Uij are even, and I{a) — 0{n~^~^) if not. 

Proof. By using Theorem 2.2 and Theorem 3.1, we obtain the following estimate: 
-^(") ^ ^^A(^i)Sa{ar)Wkn{-r^,a) 



k 



n 



n 



-k 



^5^{ai)5a{ar){5T,a + 0{n ^)) 
(#W5.(a/) = 5^(«r) = 1} + O(n-i)) 



In order to count the partitions appearing in the set on the right, let us go back to the 
multi-indices ai,ar described in the proof of Theorem 2.2. It is convenient to view both 
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these multi-indices in a rectangular way, as follows: 
(\...\ ... 1^.A\ 



p . . .p 



p . . .p 



\ 



api 



J 



q...q\ 



ail 



aiq 



p.. .p 



J 



In other words, the multi- indices ai/ar are now simply obtained from the matrix a by 
"dropping" from each entry aij a sequence of aij numbers, all equal to i/j. 

These two multi-indices, now in matrix form, have total length 2k = T^aij. We agree to 
view as well any pairing of {1, ... , 2k} in matrix form, by following the same convention. 

With this picture, the pairings vr which contribute are simply those interconnecting 
sequences of indices "dropped" from the same a^j, and this gives the following results: 

(1) In the case where one of the entries a^j is odd, there is no pairing that can contribute 
to the leading term under consideration, so we have /(a) = 0{n~''~^), and we are done. 

(2) In the case where all the entries a^j are even, the pairings that contribute to the 
leading term are those connecting points inside the pq "dropped" sets, i.e. are made out 
of a pairing of an points, a pairing of points, and so on, up to a pairing of apg points. 
Now since an x-point set has xU pairings, this gives the formula in the statement. □ 

4. Geodesic expansion method 

We have seen how to compute the asymptotic behavior of I (a), by using basic estimates 
on the Weingarten matrix. However, our results so far are effective only in the "non- 
degenerate" case, when all the entries of a are even numbers. In the general case, the 
computation of the asymptotic sign of I{a) requires a finer knowledge of the Weingarten 
matrix, notably with the exact computation of its leading term in n~^. 

In addition, the problem of computing the higher order terms, which are usually re- 
quired for delicate applications of the Weingarten formula, appears as well. 

We investigate these questions in this section and in the next one. We use a method of 
Collins [6], further processed by CoUins-Sniady p|. Let us begin with a key definition. 



Definition 4.1. The Brauer space Dk is defined as follows: 

(1) The points are the Brauer diagrams, i.e. the pairings of {1,2, 

(2) The distance function is given by d{TT, a) = k — loops{TT, a). 



2k}. 



It is indeed well-known, and elementary to check, that d satisfies the usual axioms for 
a distance function. This actually comes from some general categorical properties of the 
Brauer diagrams, which are valid in much more general situations. See [1], [9j. 

The Brauer space, which will play an important role in what follows, is by definition a 
metric space having (2A;)!! = 1.3.5 . . . {2k — 1) points. An interesting question is to find a 
"geometric" realization of this space. This will be discussed later on. 
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The series expansion of the Weingarten function in terms of paths on the Brauer space 
was originally found by Collins in [B] in the unitary case, then by Collins and Sniady [S] 
in the orthogonal case. We present below a slightly modified statement, along with a 
complete proof, by using a somewhat lighter formalism. 

Theorem 4.2. The Weingarten function Wkn has a series expansion in of the form 



9=0 



where the objects on the right are defined as follows: 

(1) A path from -n to a is a sequence p = [vr = tq 7^ ri 7^ . . . 7^ = a]. 

(2) The signature of such a path is + when r is even, and — when r is odd. 

(3) The geodesicity defect of such a path is g{p) = T,^^^d{Ti^i, Ti) — d{n, a). 

(4) Kg counts the signed paths from tt to a , with geodesicity defect g. 

Proof. Let us go back to the proof of Theorem 3.1. We can write Gkn = n^^{I + H), and 
in terms of the Brauer space distance, the formula of H is simply: 




for TT = O" 

^^^) for TT ^ a 



Consider now the set Pri^c, o") of r-paths between vr and a. According to the usual rule 
of matrix multiplication, the powers of H are given by: 

ff'-(7r,a) = Yl H{ro,n)...H{Tr-i,Tr) 

pGPr(7r,cr) 
pGPr(7r,(T) 

Thus by using the formula (1 + H)~^ = 1 — H + — + . . ., we obtain: 



r=0 

p€Pri-iT,a-) 



Now by rearranging the various terms of the double sum according to their geodesicity 
defect g = g{p), this gives the formula in the statement. □ 

For the J(a) translation of the above result, it is convenient to use the total length of 
a path, defined as d{p) = T,l^^d{Ti-i, Ti). Observe that we have d{p) = d{7r, a) + g{p). 
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Theorem 4.3. The integral I{a) has a series expansion in rT^ of the form 

oo 

J(a) =n-^^Hd{a)n-'^ 

d=0 

where the coefficient on the right can be interpreted as follows: 

(1) Starting from a G Mpxg(N), construct the multi-indices ai,ar as usual. 

(2) Call a path "a- admissible" if its endpoints satisfy 6^^(01) = 1 and 6(j{ar) = 1. 

(3) Then Hd^a) counts all a-admissihle signed paths in D^, of total length d. 

Proof. We combine first Theorem 2.2 with Theorem 4.2: 
^(o) = ^5^{ai)5„{ar)Wkn{T^,a) 

■K,<T 

00 

7r,o- g=0 

Let us denote by H^^ti, a) the number of signed paths between vr and a, of total length 
d. In terms of the new variable d = d{7T, a) + g, the above expression becomes: 

00 

/(a) = n-''Y6^{ai)6^{ar)YHd{TT,a)n-'^ 

TT,a d=0 

d=0 \ TT,a ) 

We recognize in the middle the quantity iJrf(a), and this gives the result. □ 

5. MOBIUS FUNCTION ESTIMATES 

In this section we derive some concrete consequences from the abstract results in the 
previous section. First, let us recall the following result, due to Collins and Sniady 0. 

Theorem 5.1. We have the estimate 

Wun{T^, a) = n-'=-'^(-''^)(/i(7r, a) + 0{n-^)) 

where fi is the Mobius function. 

Proof. We know from Theorem 4.2 that we have the following estimate: 

WW^, a) = n~^-'"^^'^\Ko{n, a) + 0{n-')) 

Now since one of the possible definitions of the Mobius function is that this counts the 
signed geodesic paths, we have Kq = fj,, and we are done. □ 



10 



TEODOR BANICA 



It is probably interesting to note here that in some more general frameworks, e.g. when 
On is replaced by its free version O^, the straightforward extension of Theorem 5.1 doesn't 
seem to hold in full generality. In fact, the computation of the first order term of the 
generalized Weingarten functions is an open, interesting question. See [T], f9], P]. 

Let us go back now to our integrals I (a). The analogue of the above result of Collins 
and Sniady, fully generalizing Theorem 3.2, is as follows: 

Theorem 5.2. We have the estimate 

I{a) = n-'^-"(")(/i(a) + 0{n-^)) 

where the objects on the right are as follows: 

(1) e(a) = min{d(7r,cr)|7r,cr G Dk,S^{ai) = 5<x(ar) = !}• 

(2) fi{a) counts all a-admissible signed paths in D/^, of total length e{a). 

Proof. We know from Theorem 4.3 that we have an estimate of the following type: 

J(a) =n-''-^{He{a) + 0{n-^)) 

Here, according to the various notations in the previous section, e G N is the smallest 
total length of an a-admissible path, and He{a) counts all signed a-admissible paths of 
total length e. Now since the smallest total length of such a path is of course attained 
when the path is just a segment, we have e = e(a) and Hf.{a) = fi{a), and we are done. □ 

Summarizing, we have now a full description of the asymptotic behavior of /(a). To 
illustrate how Theorem 5.2 applies, let us recover the results in Theorem 3.2: 

(1) Assume first that all entries of a are even. In this case there is at least one partition 
TT such that St^^oi) = 5^(ar) = 1, so we have e(a) = 0, and we recover the leading 
term from Theorem 3.2. As for the coefficient /i(a), this counts the partitions 
71 such that 5,r(oz) = Sn{(ir) = 1, SO wc fuUy recover Theorem 3.2. 

(2) Assume now that a has at least one odd entry. Here there is no vr as above, so we 
have e(a) > 1, and we recover the quantity 0{n~''~^) from Theorem 3.2. 

Regarding now the higher order terms, the situation is quite complicated. In fact. 
Theorem 4.3 is quite difficult to use as stated, as it was the case in fact with Theorem 
4.2 as well. Let us point out, however, that a positive answer to the following question 
would definitely change the situation, by making these two results fully efficient: 

Question 5.3. Is there a "geometric" model for the Brauer space? 

In other words, what we are asking for is a better model, of geometric or perhaps 
analytic nature, which would convert the quite heavy counting of geodesies on the Brauer 
space into a kind of calculus computation. As far as we know, there is no idea yet here. 
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6. Zonal function approach 

We have now a quite satisfactory picture of lim^^oo I {a)- However, the computation of 
/(a) for fixed n remains a quite subtle problem, as shown for instance by Theorem 1.3. 

As explained in the previous section, the exact formula coming from the geodesic ex- 
pansion is quite difficult to handle, in the lack of a better understanding of the Brauer 
space Dk- However, there have been many attempts for working out the combinatorics of 
the geodesic expansion. Let us mention here the pioneering work of Collins in the unitary 
case [6], further processed by Collins and Sniady in the orthogonal case [8]. 

A remarkable advance on this subject came quite recently, in the preprint of Collins 
and Matsumoto f7] . In this section we briefly explain their new formula, and we work out 
the corresponding consequences regarding the integrals /(a). 

The formula of Collins and Matsumoto [7] (see also Zinn- Justin [13]) is as follows. 

Theorem 6.1. We have the formula 

'^"'"•"' = (2*)"n,„,«(« + 2,-.-i) 

where the various objects on the right are as follows: 

(1) The sum is over all partitions o/{l, . . . , 2k} of length /(A) < k. 

(2) is the corresponding zonal spherical function of {S2k, Hk) ■ 

(3) ^■s the character of S2k associated to 2A = (2Ai,2A2, • . .). 

(4) The product is over all squares of the Young diagram of X. 

It is of course possible to deduce from this a new a formula for /(a), just by putting 
together the formulae in Theorem 2.2 and Theorem 6.1. However, there are probably 
several non-trivial simplifications that might appear when doing the sum over vr, a, and 
we do not know how the final statement about I{a) should look like. 

Instead, let us just record the following consequence. 

Corollary 6.2. The possible poles of I {a) can be at the numbers 

-{k - 1), -{k - 2), ... ,2k - 1,2k 
where A; G N associated to the admissible matrix a G Mpxg(N) is given by k = Sajj/2. 

Proof. We know from Theorem 2.2 that the possible poles of I{a) can only come from 
those of the Weingarten function. On the other hand. Theorem 6.1 tells us that these 
latter poles are located at the numbers of the form —2j + i + 1, with ranging over 
all possible squares of all possible Young diagrams, and this gives the result. □ 

As explained in [7], the main feature of the zonal function approach is its remarkable 
efficiency in respect with numeric computations. It is our hope that the above results 
would be of help in producing, perhaps via a computer program, more "data" for the 
exact computation of /(a). See the comments in the concluding section below. 
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7. Concluding remarks 

We have seen in this paper that the representation theory approach to the computation 
of polynomial integrals over On, originating from Weingarten's paper [12], formulated 
and developed by Collins and Sniady in [B], [S], further processed in [7|, [TT], [13], and 
helped with some new ideas from [T], [5], [ID], leads to a number of concrete results on 
the integrals of type /(a), which can be regarded as being the main objects of study. 

Together with some explicit new computations, to be presented in our forthcoming 
paper [3], and based this time on purely geometric methods, these results give hope for 
the existence of a general explicit correspondence of type a /(a), with the output 
/(a) G Q{n) depending somehow "arithmetically" on the input a G Mpxg(N). 

We do not know if it is so, but we intend to intensively pursue the exploration of this 
correspondence, which would be of course extremely useful for any kind of mathematical 
area where the Weingarten formula has already proved to be useful. Note in particular 
that the present results provide a much better hope for the program launched in [2j. 
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